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Abstract: In this paper, we present an exactly solvable phase transition model in which 
the phase transition is purely statistically derived. The phase transition in this model is 
a generalized Bose-Einstein condensation. The exact expression of the thermodynamic 
quantity which can simultaneously describe both gas phase and condensed phase is solved 
with the help of the homogeneous Riemann-Hilbert problem, so one can judge whether 
there exists a phase transition and determine the phase transition point mathematically 
rigorously. A generalized statistics in which the maximum occupation numbers of different 
quantum states can take on different values is introduced, as a generalization of Bose- 
Einstein and Fermi-Dirac statistics. 
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1 Introduction 

A few exactly solvable models play important roles in phase transition theory, since most, 
if not all, of our understanding of phase transitions comes from studying models In 
this paper, we present a purely statistically derived solvable phase transition model. In the 
model, the exactly solved thermodynamic quantity can simultaneously describe different 
phases. Therefore, whether there is a phase transition can be judged mathematically 
rigorously, and the phase transition temperature can be calculated directly by analyzing 
the discontinuity in the thermodynamic quantities or their derivatives. 

Bose-Einstein condensation (BEC) is the first purely statistically derived example of a 
phase transition. The phase transition in the present model is a generalized Bose-Einstein 
condensation; in other words, the phase transition is a BEC type phase transition. 

The BEC type phase transition is a sudden change in the microscopic particle dis- 
tribution: in the gas phase, no quantum state is macroscopically occupied, while in the 
condensed phase, there is a quantum state being occupied by a macroscopic number of 
particles. The microscopic particle distribution determines the macroscopic behavior of a 
thermodynamic system, or, the macroscopic behavior reflects the average contribution of 
all quantum states in the system. In the condensed phase, the macroscopic behavior of 
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the system is to a certain extent determined by the single quantum state that is macro- 
scopically occupied, since the number of particles in such a state is of the same order of 
magnitude of the total number of particles. As a macroscopic manifestation of such a 
sudden change in the particle distribution, there is a singularity in the thermodynamic 
quantity. 

We will show that, beyond the Bose case, there are still other systems that can display 
BEC type phase transitions. 

First, we will introduce a generalized statistics in which different quantum states have 
different maximum occupation numbers, and Bose-Einstein and Fermi-Dirac statistics are 
its special cases. Especially, we will pay attention to a special case of the generalized 
statistics in which at least one quantum state's maximum occupation number is infinite 
and show that the BEC type phase transition may occur in such systems. For example, 
we will show that a BEC type phase transition can occur in an any-dimensional ideal gas 
obeying the generalized statistics in which the maximum occupation number of the ground 
state is infinity, like that in the Bose-Einstein case, and the maximum occupation number 
of all other quantum states is 1, like that in the Fermi-Dirac case. For comparison, recall 
that the BEC can occur only in three-dimensional ideal Bose gases, but cannot occur in 
one- and two-dimensional cases. 

The mathematical method for solving the model is based on the homogeneous Riemann- 
Hilbert problem — the boundary problem of analytic functions, which comes from the 
theory of singular integral equations [21 [3] . 

Moreover, our result also shows that a phase transition occurs only in the thermo- 
dynamic limit, i.e., the total number of particles and the volume must be infinite. The 
common proof for this result depends on an assumption that a finite volume can accom- 
modate at most a finite number of particles, which is, of course, only valid for non-ideal 
gases jl]. Our result provides an example that this result holds also for ideal gases. 

In section 2, we introduce the generalized statistics. In section 3, we construct and 
solve the phase transition model. Discussions and an outlook are given in section 4. 



2 The generalized statistics 

In this section, we introduce a generalized statistics in which the maximum occupation 
number of a quantum state can take on unrestricted integers or infinity and the maximum 
occupation numbers of different states may be different. 

Let TLi be the maximum occupation number of the i-th quantum state, where tii can 
take on an integer or oo. The grand partition function is 

~ 1 _ e-/3K+i){e-/^) 

1=0 

where T is the temperature, V the volume, /i the chemical potential, Ei the energy of the 
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i-th state, and (3 = 1/ (kBT). Then the equation of state reads 
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The equations of state for Bose-Einstein, Fermi-Dirac, and Gentile [SJ E] cases can be 
recovered by setting m = oo, rii = 1, and rii = n, respectively. 

In this paper, we consider an ideal gas obeying the generalized statistics in which the 
maximum occupation number of only one quantum state is oo, but of all other quantum 
states is n, a given integer, i.e., = oo and = n, (i ^ k). The equation of state for such 
an ideal gas with the dispersion relation e = / (2m) in i/ dimensions can be obtained 
from equations ([2]) and 
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where z = e^^ is the fugacity, Nx = ^p||^y^j ^ , and A 



[2n^/2mkBT) 



17^ 



(4) 
(5) 

is the mean thermal 



wavelength. ha{z) can be expressed by the Bose-Einstein integral, ga{z), as 

K{z)=g,{z)-{n+l)-^''-^^g,{z^^+^), 

and in the limit n ^ 00 or n = 1, h„{z) returns to the Bose-Einstein integral Qaiz) or the 
Fermi-Dirac integral faiz), respectively 0. 

We will show that such an ideal gas system may display the BEG type phase transition, 
and whether the phase transition occurs or not lies on the value of k, the position of the 
state with an infinite maximum occupation number in the spectrum. 



3 The phase transition 

In this section, we consider two cases which can display BEG type phase transitions: the 
ideal gases obeying the generalized statistics with tiq = 00, rii = n {i ^ 0) and with 
Hk = 00 {k ^ 0), Ui = n (i ^ k). An interesting case is n = 1. In this case, the maximum 
occupation number of only one state is the same as that in the Bose case, but of all other 
states is the same as that in the Fermi case. We will show that even such systems in which 
only one state's maximum occupation number is infinite can still display BEG type phase 
transitions. 
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3.1 The explicit expression for the fugacity 

For judging whether there is a phase transition or not and determining the phase transition 
temperature, we first solve the exact exphcit expression for the fugacity from the equation 
of state, and, then, analyze the discontinuity in the derivative of the fugacity. 

Based on the homogeneous Riemann-Hilbert problem [3], we can solve the explicit 
expression for the fugacity from equation ([5]) exactly (A brief introduction to the method 
of the Riemann-Hilbert problem see Ref. [7]). 

For the case of = oo and rii = n (i ^ k), introduce a complex function 

where (C) is an analytic continuation of h„{z). From equation ([5]), we can see that the 
fugacity is a zero of ^ {C) on the real axis. Therefore, the problem of solving the fugacity 
z is converted into the problem of seeking the real zero of the complex function ^ (Q. 
We can express ^' (C) as 

(c-^) n (^-^^) 
n (c - n (c - cmT- 

j=l m=l 

where (p (C) is the fundamental solution of the homogeneous Riemann-Hilbert problem, 
which has no zeros and singularities, z (the fugacity) and are zeros of ^ {(), Uzero is 
the number of the zeros, pj is a pole of ^' (C), ^^p is the number of the poles, Cm is an 
endpoint that is different from infinity of the boundary of the analytic region of ^ {Q (in 
the present case, the boundary of the analytic region of ^ ((") is a set of rays (see figure 
1) with the origins Cm)-, n^, is the number of the endpoints different from infinity of the 
boundary (in the present case is the number of the rays), the constant Km is introduced 
to equal the degrees of divergence of the two sides of this equation at the m-th endpoint 
Cm, and r/ is a constant. 

From equation ([7|) , it is easy to see that we can in principle obtain an explicit expression 
of the fugacity z. For this purpose, we need to first determine the fundamental solution 
99(C), the endpoints Cm, and the poles pj, etc. 

The analytic region. For determining Cm , the endpoints of the boundary of the analytic 
region, we first analyze the analytic region of ^ {Q). 

The boundary of the analytic region of ^ {Q is determined by the analytic region of 

K (C) = 9. (C) - - — \-^,9. {C^') , (8) 

[n + 1) 

where {Q is the analytically continued Bose-Einstein integral which is just the Jonquiere 
function [8]: 

g. (C) = Li^ (C) • 
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The boundary of the analytic region of Liu (C) is the positive real axis from 1 to oo [8]. 
Consequently, the boundary of the analytic region of (C) and ^ {() consists of n rays 
with origins on the unit circle (figure 1), denoted as Lfn, m = 1,2,- •• ,n, i.e., the m-th 
ray is 



■ 27rm 



e*"+i , ooe "+i j . It should be emphasized that (C) has no singularity on the 
positive real axis, or, ((") is analytic on the positive real axis. Then, the endpoints that 
are different from infinity of the boundary (the origins of the rays) are 



■ 27rm 



Cm = e*"+i, m = l,2,---,n. (9) 

Therefore we have rib = n. 

The fundamental solution of the homogeneous Riemann-Hilbert problem. Now, we 
calculate the fundamental solution of the homogeneous Riemann-Hilbert problem, <^(C)- 

Using the result of the homogeneous Riemann-Hilbert problem [3] , we have 

n 

<^(C) = e^(^)n (C-Cn^)''". (10) 

771=1 

Here 

InG(x) 



dx- 



x-C 



■ 2-Km 

" , roo InG xe "+i 



where 



1 ^ ^ -27™ 1,1/ , , 

^ E «^ / ^ (11) 

G(C) = ^ (12) 

is the jump of (C) on the boundary, and Am is an integer determined by the condition 

InG(cm) , , „ InG(cm) . . , 
^ Re h Am = U, II =F Re is an integer, 

-1< T Re ^"'^(^"^) + Am < 0, otherwise. (13) 
27ri 

We first need to analyze the analytic region of the fundamental solution <^(C)- From 
equation ([7]), we can see that the boundary of the analytic region of (C) consists of the 
non-isolated singularities of ^ (C), and the jump of (C) on the boundary is the same as 
that of ^' (C), i.e., 

G(C) = ^. (14) 

The value of ^ (C) on the two sides of Lm is determined by the value of the function 
(C) on the two sides of L^: 



/i^ xe "+i = hce "+i - -— ^g^ (x"+ ) =F ^tttt (l'^^;) , m = l,2,--- ,n, 

\ / \ / (^72 -|- 1) 

(15) 
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where (C) is the Cauchy principal value of {() on the boundary [9]. Then 
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1 =F *TT7 — rT~ (Inx) ' 



r(i//s) 



(16) 



Note that, in the case of = oo and rii = 1 (i ^ k), equation (I16p reduces to 



= (-X) + 



(17) 



where f,^/^ (^) is the Cauchy principal value of the analytically continued Fermi-Dirac 
integral. 

Next, we calculate Am from equation (flSll . 

At the endpoints of the boundary (including both the endpoints that are different from 
infinity and the infinity), we have 



G[e — ) =G{oo) = 1. 



Choosing InG (oo) = gives 



In G e '^ = i arg G e ~ 



-i27r. 



Then we have 



Therefore, the fundamental solution is 



=7(0. 



C-1 

Qn+l _ ■ 



(18) 

(19) 
(20) 

(21) 



The value of Km- The parameter is chosen to guarantee the degrees of divergence 
of the two sides of equation d?]) at the endpoint Cm to be the same. 

At the origin of L^, Cm = e*2™/("-+i)^ when v/s > 1, g^/g (x"+^) and (Inx)"^''^^ are 
convergent, and when u/s < 1, the degrees of divergence of g^^/g (x""*"^) and (Inx)'^''''"^ 
are less than one. Thus we have 

Km = 0. (22) 

The isolated singularity. ^ ((") has only one isolated singularity, 

p = e'^'K (23) 
The number of zeros of ^' {(). Substituting the above result into equation d?]) gives 



ra+l 



1 



C-1 



(24) 
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Figure 1: The analytic region of ^' (C)- The cases of n = even and n = odd are illustrated 
in (a) and (b), respectively. 

In principle, if and rj are known, one can obtain the explicit expression of z directly. 
Nevertheless, the difficulty of finding the zeros is often the same as the difficulty of 
finding the zero z. Alternatively, we can construct a set of equations for z, uJi, and rj, and 
obtain z by solving the equations. 
For solving z, we need 

'^zero~^^ equations. Since the number of the isolated singularities 
of ^ ((^) is already known, the number of the zeros, Uzero, can be determined by the 
argument principle, the contour being chosen as in figure 1. The result shows that ^ (^) 
has n + 1 zeros on the C-plane, so we need n + 2 equations for determining z. 

The case of rik = oo and Ui = 1 (i ^ k). For simplicity, we consider the case of 
Uk = oo and rii = 1 (i ^ k); the solutions for more general cases can also be obtained 
exactly but with more complex forms. When n = 1, equation (j24p becomes 



substituting equation ([26]) into the derivative of equation (|25]) and setting C = give 



rj iC - z){C- CO) = e-''(f) (C + 1) (C - 6^^*=) * (C) 



(25) 



In this case, for solving z, we need three equations. 
Substituting C = into equation ((25]) gives 

7]ZL0 = e^'^-"'^°^; 



(26) 



y (0) - 1 + e-^^^ - i - - 



(27) 



Z OJ 

substituting ( = e^'^'' into equation ([25]) gives 



N 




(28) 



8 



Solving the equations (j26|) . (j27[) and (j28|) . we have 



(29) 



where r?A = iVA/iV + 7' (0) - 1 and 



1 , InG(-x) , , 

3.2 The phase transition and the necessary condition for phase transi- 
tions — the thermodynamic Hmit 

Equation (|29p is an exact expression of the fugacity, which simultaneously describes both 
gas phase and condensed phase, and, of course, can describe the transition between these 
two phases. From equation (I29p . we can directly see that the thermodynamic limit, N — > 
00, is the necessary condition for the phase transition. 

The fugacity given by equation (j29p is a smooth function and there is no singularity. 
Therefore, there is no phase transition regardless of how low the temperature is. However, 
in the thermodynamic limit, i.e., N ^ 00, equation ()29p becomes 

1 J 6^^*=, when rjx < e~f^^*' , 



r]\ + e ^^fe + \r]\ - e ^^*| ^, when rjx > e . 



(31) 



The discontinuity may appear in the first-order derivative of the fugacity and the phase 
transition may occur. The discontinuous point appears at 

Vx = e-^'\ (32) 

which is just the phase transition point. 

In the phase transition theory, there is a fundamental law: the necessary condition for 
a phase transition is that the system must be infinite, i.e., the thermodynamic limit. The 
proof of this statement depends on the assumption that a finite volume can accommodate 
at most a finite number of particles. If the number of particles is finite, the partition 
function will be an analytic function and, consequently, there is no singularity in the 
thermodynamic function and there is no phase transition ^[TnilHI- Such an assumption 
is equivalent to requiring that the particle must have a nonzero volume. Clearly, this 
assumption does not hold for ideal gases. That is to say, though this conclusion is valid 
for all realistic systems (realistic gases are non- ideal gases), this proof is not valid for the 
idealized model — ideal gases. Our above result shows that for ideal gas systems, the 
thermodynamic limit is still a condition for a phase transition. 

3.3 The case of tiq = C)0 and rii = 1 (i ^ 0): the phase transition temper- 
ature 

We first consider the case of no = 00 and rii = 1 {i ^ 0), i.e., the state whose maximum 
occupation number is infinite is the ground state, Sk = £0 = 0. 
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In any dimension, there must exist a phase transition. This can be verified directly by 
observing the discontinuity in the derivative of the fugacity z from equation (|32]) and the 
transition temperature is determined by 

r/A = 1. (33) 
Then the phase transition temperature reads 



N sT {u/2) 1 1 

V2T{u/s){l-2^~^/')C{i^/s)\ ' ^ ' 



Now let us see what happens when a phase transition occurs. The total number of the 
excited particles, from equation (0), is 

Ne = N^fu/s (z) ■ (35) 

Comparing the expressions of Tc and gives that when the phase transition occurs, 

Ne = N. (36) 

This is just the condition that one determines the phase transition temperature for a 
BEC in an ideal Bose gas. In our case, however, this result comes from a mathematically 
rigorous calculation rather than being put in by hand. 

This result indicates that when the phase transition occurs, the macroscopic properties 
of the system will begin to be controlled, to a certain extent, by a unique quantum state 
(here the state is the ground state). Such a phase transition is a sudden change in the 
particle distribution: in the gas phase, the macroscopic behavior of the system is a mean 
contribution of all quantum states, but in the condensed phase, the quantum state with 
an infinite maximum occupation number dominates. This is a BEC type phase transi- 
tion. More concretely, when the phase transition begins, the number of excited particles 
decreases as the temperature decreases, while the number of particles in the ground state 
increases as the temperature decreases: 

n-[tJ ' iv - ^ [tJ ■ ^^^^ 

Different from the BEC in an ideal Bose gas, the BEC type phase transition can occur 
in any dimension in the ideal gases obeying the statistics in which the maximum occupation 
number of the ground state is oo and of all other states is finite, since the phase transition 
can occur for any positive value of u/s. In the Bose case, however, the BEC only occurs 
when v/s > 1, and, as a result, the BEC only occurs in three-dimensional Bose gases. 
This is because the Bose-Einstein integral in the Bose case is replaced by /i^/g (z) in the 
present case, while h^/g (z) is always bounded for < z < 1. That is to say, in such an 
ideal generalized-statistics gas, the occurrence of the phase transition is easier than that 
in a Bose system. 
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3.4 The case of rik = oo and rij = 1 {i ^ k): the phase transition temper- 
ature and the Fermi energy 

Next, we consider the case that the only state with an infinite maximum occupation 
number is not the ground state, i.e., = oo and rii = 1 {k ^ and i ^ k). More general 
cases can be treated by the same procedure. 

The phase transition temperature. Prom equation (|31|) . we can see that the discon- 
tinuous point of the derivative of the fugacity appears at rjx = e~^^*=. Equation (f3TI) 
indicates that the phase transition appears at z = w = 6^^*=. Substituting ( = e^^'^ into 
the derivative of equation ([25]). when N ^ oo, gives 

A(fcijT,r/V./.(e^) =1, (38) 

where A = ar(^y2j h^'^ — Tf' This result indicates that when the phase transition 
occurs, the number of particles in all the states except the infinite-maximum-occupation- 
number k-th. state equals the total number of particles of the system, i.e., A^kt^oo = 

Based on the homogeneous Riemann-Hilbert problem, we can solve the phase transition 
temperature from equation ()38p . For simplicity, we only give the result for the case of 
s = v. 

Introduce a complex function 

-kBTh e^B- -1. (39) 



' ' sr(s/2)iV 

The phase transition temperature is a zero of (r) on the real axis. 

The analytic region. We first analyze the analytic region of (p (r) on the r-plane. The 
analytic region of (p (r) is determined by the behavior of the analytically continued Fermi- 
Dirac integral, /i (e^'=/(^sT)) ^ which is illustrated in figure 2(a). The boundary of this 
region is complex. Introducing a transformation 

^ = (40) 
ksT 

we have 

^(6 = ^^-^-/ife-^l-l (41) 
sTis/2)N ^^'\ J ^ ' 



The boundary of the analytic region of ^ (^) on the .^-plane is 
ReC > 0, 

Ime = ^^^±^, , = 0,±1,±2,..., (42) 

as illustrated in figure 2(b), which is a set of rays running parallel to the real axis with 
origins 

(2fl+l)7r\ , , 

0, ^ ^ . (43) 
J 
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Figure 2: (a) The analytic region of (P{t); (b) The analytic region of "0 (C)- 

The fundamental solution of the homogeneous Riemann-Hilbert problem. We can also 
express ip {$,) in the form of equation ([7]), and, then, solve the explicit expression of the 
phase transition temperature. First, we seek for the fundamental solution of the homo- 
geneous Riemann-Hilbert problem. According to equation (fTOj). the fundamental solution 
can be written in the following form: 

oo 

<^(^) = e^«) n iC-Cgf", (44) 

g=— oo 

where 

7K) = 5^-/ ^x'^, (45) 
and the integral is along the boundary of the analytic region, 

Lg:^ = x + i (^i±l)Z^ X G [0, oo) and q = 0, ±1, ±2, • • • . (46) 

The constant is an integer satisfying the condition p3|) . 

The jump on the boundary of the fundamental solution (/? (^) is the same as that of 

CO^^^^lf^ ,47) 

f (4) (0 
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■0^ (0' value of (^) on the two sides of the boundary, is determined by the behavior 
of the analyticahy continued Fermi-Dirac integral, 



(48) 



Then, 



■0 \ X + i 



(2g + l)7r 



A' 



(2g+l)7r 



(2g + l)^ 



A' 



(2g+l)7r 



(2g+l)7r 



(49) 



where A' 



2 27r«/2m V 



^P^^^2) — p — TV' From equations (f^7|l and (H9]l . we can see that G (oo + z (2g + 1) vr/e/c) 
1. The constant \q is determined by the condition (fT3]) . Choosing In G(cxd) = gives 



=0. 



Consequently, the fundamental solution is 



(50) 



(51) 



The value of Kg. At the endpoints Cg = (0, (2g + 1) vr/cfc), we have -0 (2(? + 1) tt/sa,.) 
ln[^-i(2^ + l)7r/efc]. Then, 



Kg = 0. 



(52) 



The isolated singularity of ip {(,)■ ip (^) has only one isolated singularity, 

p = 0. (53) 

The number of the zeros of ^ (^). By the argument principle, the contour being illus- 
trated in figure 2(b), we can determine that V (0 has only one zero, ^ = f3c = 1/ {kBTc), 
which is on the real axis. 

Introducing ^ (C) = ■^'0 (C)> '^^ have 



<I>(e) = i;(C-/3c)e^«), 



(54) 



where f is a constant. Substituting ^ = into equation ()54p and its first-order derivative 
gives two equations. Solving these equations gives 



1 



1 efc 



2Tr' /'^mkB V 2 V27 ln2 21n2A;ij 
1 1 
ks 2Tii 



g=— oo 



°°^^lnG(x + i(2g + l)7r/efc) 



{x + i{2q + l)Ti/ekY 



(55) 



The last term of equation (|55p is small when is small, and is roughly proportional to 
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The explicit expression of the phase transition temperature shows that depends on 
the value of Ek- Especially, Tc = appears at 



i.e., when > ep, there will be no phase transition. Note that this result holds also for 
the case oi s ^ v. It is not difficult to recognize the physical meaning ol ep'- it is just the 
Fermi energy of a z/-dimensional ideal Fermi gas with the dispersion relation e = p**/ (2m) 
[5]. The reason why there is no phase transition when > ei;' is that if > ep^ the 
states below ep can accommodate all particles in the system and, then, there are no enough 
particles accumulating in the fe-th state, i.e., the BEC type phase transition cannot occur. 

4 Discussions and outlook 

In this paper, we construct an exactly solvable phase transition model. We first consider 
a generalized statistics in which the maximum occupation numbers of different quantum 
states can take on different values. When the maximum occupation numbers of all the 
states are the same, e.g., equaling oo, 1, or an arbitrary integer, the generalized statistics 
returns to Bose-Einstein, Fermi-Dirac, or Gentile statistics [5l [6l [12], respectively. The 
model constructed in this paper is an ideal gas obeying the generalized statistics in which 
the maximum occupation number of only one state is infinite, but of all other states is 
finite. The phase transition which occurs in such systems is the BEC type phase transition. 
For judging if the phase transition can occur and determining the phase transition point, 
we calculate the exact explicit solution for the fugacity with the help of the mathematical 
result of the homogeneous Riemann-Hilbert problem. By observing the discontinuity in 
the derivative of the fugacity, we analyze the phase transition rigorously. From this phase 
transition model, we can see that the thermodynamic limit is a necessary condition for a 
phase transition of an ideal system. 

For constructing the solvable phase transition model, we introduce a kind of interme- 
diate statistics. Various generalized exclusion statistics models play important roles in 
many fields [Gj [13l [HI [151 [13 [13 [18], since many physical systems may behave as neither 
Bose-Einstein nor Fermi-Dirac system. Though nature realizes only bosons and fermions, 
there are many composite-particle systems, e.g., the Cooper pair in the theory of super- 
conductivity, the Fermi gas superfluid [19], the exciton [20], the magnon [21], etc. For 
example, a boson consists of two fermions obeys Bose-Einstein statistics. However, when 
two such bosons come closer together, the fermions in the composite bosons may "feel" 
each other, and the statistics may somewhat deviate from Bose-Einstein statistics. In this 
case, such a composite system can be effectively viewed as obeying a kind of intermediate 
statistics. It is shown in a resent study [22] that the fermion pairs in the one-dimensional 
Fermi gases obey generalized exclusion statistics. 




(56) 
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The experimental and theoretical research of BEC is a branch of the statistical physics 
of a rapidly growing importance [23]. The BEC of ideal Bose gases is a special case of 
the generalized BEC phase transition. By studying this exactly solvable model, we can 
also obtain a deeper insight into the BEC of ideal Bose gases. We can conclude that the 
conditions for the BEC type phase transition as follows: 

(1) There must exist a low enough quantum state with an infinite maximum occupation 
number, where "low enough" means that when the temperature tends to the absolute 
zero, there are still a macroscopic number of particles in this state. In other words, this 
condition requires that the total capacity of all states below such a state must be small 
enough so that this state can be macroscopically occupied when the temperature is low. 
As a result, the energy of this state, denoted as e™™, must be the lowest one among the 
states whose maximum occupation numbers are infinite. Such a condition is of course 
satisfied by a Bose system since the maximum occupation number of the ground state is 
infinite. However, for the systems obeying the generalized statistics, as discussed above, 
in the case of no = co and rii = n (i 7^ 0), this condition can be always satisfied, but in 
more general cases, e.g., the case of = 00 and rii = n {i ^ k), this condition can be 
satisfied only when < ep. 

(2) The state that will be macroscopically occupied when a BEC type phase transition 
occurs must be isolated from other infinite-maximum-occupation-number states, where 
"isolated" means that the state density of this state is a (^-function, i.e., the state density 
of the states with infinite maximum occupation numbers (not the state density of the 
system) must take the form of poo (e) = Voo (e) + 6 [e — and r/00 (^S)™) = 0; where 
Poo (e) is the density of the infinite-maximum-occupation-number states. In the examples 
of the generalized statistics we considered above, this condition is satisfied naturally, since 
there is only one infinite-maximum-occupation-number state, for no = 00 and Ui = n 
(i 0), Poo (s) = (^(e), and for Uk = 00 and rn = n {i ^ k), poo (e) = S{e-ek)- 
However, for the case of ideal Bose gases, this condition is not always satisfied. In ideal 
Bose gas systems, the maximum occupation number of all states is infinite, i.e., the state 
density of the system p{e) = poo{£), and then the lowest infinite-maximum-occupation- 
number state is the ground state, i.e., e™™ = 0. In three dimensions, the state density is 
Poo (e) = i]oo (e) + (5 (s)) where r/00 (e) oc ^/e, so r/00 (0) = 0. The condition is satisfied, and 
the BEC phase transition can occur in a three-dimensional ideal Bose gas. In one and two 
dimensions, the state densities are poo (e) = 1I00 (e) oc l/^/e and poo (e) = "Hoc (s) = const, 
respectively; the above condition is not satisfied, and there are no BEC phase transitions 
in one- and two-dimensional Bose gases. 

Furthermore, many physical systems possess other kinds of statistics beyond Bose- 
Einstein and Fermi-Dirac statistics. For example, the Calogero-Sutherland model is shown 
to possess fractional statistics [15) . a spinless fermions system in two dimensions may 
obey exclusion statistics |17) . and bound pairs of fermions form hard-core bosons obeying 
generalized exclusion statistics [22]. Moreover, in the model constructed in the present 



15 



paper, there are both bosonic and fermionic states in a system. In a Bose system, if 
each boson consists of two fermions, then in the system there must simultaneously exist 
both bosons and fermions due to the fact that there exists an "ionization" energy. As 
long as the temperature of the system is not the absolute zero, there are always a certain 
proportion of particles having energies larger than the "ionization" energy and behaving 
as fermions. In such a case, the particle at the low-lying state behaves as a boson and the 
particle at the high-lying state behaves as a fermion. That is to say, a composite system 
will not accurately possess Bose-Einstein or Fermi-Dirac statistics. In such a composite 
system, our model may work. We will address this issue in future work. Moreover, a 
system consisted of both bosons and fermions also has been studied in literature |24j . 
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